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Abstract 

In this work we extend simple top-hat model of structure formation to the two- 
component system made of baryonic and dark matter. We use Harrison-Zeldovich 
spectrum as the initial condition for the structures and calculate their evolution up 
to the present time. While we do not take into account some complications during 
the structure formation, such as the merging of galaxies, however this formalism 
can give us a qualitative picture from the formation of structures. We show that in 
this model small scale structures evolve faster than the larger ones and it predicts a 
down-top scenario for the structure formation. The trend of power spectrum in this 
model is compatible with the observations and results in ~ 0.8. This formalism 
provides an analytic treatment of structure growth and can easily show the effect 
of the cosmological parameters on the formation of the structures. As an example, 
the effect of a parameterized dark energy model on the growth of the structures is 
investigated. 
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1 Introduction 



Observations of tlie Cosmic Microwave Background (CMB) by the COBE 
satellite and the subsequent experiments such as WMAP indicate the ex- 
istence of temperature fluctuations of order ~ 10~^ at the last scattering 
surface (El; Q; y)- This temperature contrast on CMB may result from the 
primordial quantum fluctuations at the early universe. The quantum fluctua- 
tions in the inflationary scenario, provide a specific spectrum for the matter 
so-called Harrison-Zeldovich. However, Recent observations by WMAP show 
a small deviation from this spectrum (us = 0.958 ±0.016) (2). An outstanding 
characteristic of this spectrum is its scale-independent property, means that 
all the perturbations have some density contrast of about 10~^ at the entering 
time to the horizon (i5i). One of the questions in cosmology is that of how these 
small perturbations at z^ec = 1100 can grow to the present non-linear struc- 
tures while we expect from the standard structure formation theory that they 
should evolve to a density contrast of 5 ~ 6dec x Zdec = 10~^ at the present 
time (|6|). 

Including dark matter as one of the components of the universe is a solution to 
this question. Since the entering time of the structures to the horizon depends 
on the size of structure, we expect small scale structures enter earlier than the 
larger ones. Meanwhile, before the decoupling epoch, pressure of the radiation 
prevents formation of baryonic structures smaller than the Silk mass (7), the 
dark matter structures continue their growth. After decoupling, the mutual 
interaction of dark matter-baryonic matter speeds up the growing rate of the 
structure formation and results in non-linear structures at the present time. 

Here in this work we use two components of baryonic and dark matter in the 
top-hat model for studying their growth. We take the initial condition of sub- 
horizon over-dense regions at the last scattering surface. The size of pertur- 
bations in our concern guarantees using of Newtonian mechanics (jsj; [q]). After 
decoupling, the top-hat sphere starts to expand up to a maximum radius and 
then turns-around to collapse. During the collapse, once the structure satisfies 
the virial condition, the global velocity turns into the dispersion velocity and 
thermalizes the gas of structure and finally prevents a catastrophic collapse 
of the structure. The result of thermalization is the ionization of the bary- 
onic gas and consequently gas starts to cool through the radiation. Cooling 
makes baryonic structure to contract further and finally the baryonic compo- 
nent reaches to a stable stage. We should point out that in this scenario we 
ignore merging effects during the formation of the structures and this model 
can be applicable only for the isolated systems. Finally we compare ag from 
this simple theoretical model with that from weak-lensing. 

The organization of paper is as follows: In Section [2] we introduce Standard 
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Collapse Model (SCM), extend it for the two component fluid and obtain 
dynamics and power spectrum of the structures. Section [3] discusses about 
the effect of variable dark energy on the formation of structure in top-hat 
model. In Section H] we study the cooling effect on the evolution of baryonic 
matter after the thermalization and estimate corresponding redshift for the 
star formation. We conclude in section [51 



2 Spherical Top-Hat Model: Structure Formation 



In this section we review the standard spherical collapse model. Here we take 
a spatial uniform distribution of the matter inside the structure, so-called top- 
hat distribution which is slightly denser than the uniform background density. 
One of the advantages of this model is that it has an analytical solution for 
the dynamics of the structure. In this section flrst we introduce the standard 
top-hat model and then extend it, introducing the two component fluid in 
the structure: (i) a non-dissipative dark matter and (ii) a dissipative baryonic 
matter. 



2.1 standard top-hat model 



A simple approach for studying the structure formation in the universe is the 
spherical collapse model. We take a spherical region embedded in the uniform 
background which has a tiny density deviation from that of the background. 
The scale of this region is much smaller than horizon and the velocities are non- 
relativistic. These two conditions guarantee the application of the Newtonian 
gravity for studying the growth of the structures (0; 0). 

For a spherical region with the radius R(t) and uniformly distributed mass of 
M, containing non-relativistic matter, the density contrast is given by: 



where phit) is the homogenous background density of the universe. The energy 
and momentum equations for a non-dissipative spherical matter is given by: 



£ = (2) 

, GM 

m=-^ (3) 
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where E is energy per mass and can be calculated from the initial condition of 
the structure, (i.e. E = Ei). The initial radial velocity of the structure is taken 
by Vi = HiRi + v'f'^''\ where Hi is the Hubble parameter of the background 
and Ri is the size of the structure at the initial time, v^^^'^^ is the peculiar 
velocity of the structure and can be given by f = —HiRiSi (@). Using the 
dependence of the peculiar velocity to the density contrast, the radial velocity 
of the structure at the initial time is: 

Vi = H,R,{l-6,), (4) 



where 5j depends on the size of the structure. Using equation (j4]), the Kinetic 
energy per unit mass of the structure is Ki = kI^\i — 25,) where Ki^"* is the 
Kinetic energy of background at a distance Ri from the center of coordinatJ^. 
For the initial potential energy of structure we have Ui = Q,iK^^\l + 6i). The 
total energy is given by the sum of the kinetic and the potential energy of the 
structure at the initial time as: 

E = -A;f a(l + 6i + 25ifiri _ ar^). (5) 

For the case of spatially flat universe, E = —?)kf^5i. Integrating from equation 
([2]) results in the equation of motion in the parametric form: 

R{e) = A{l-cos9), (6) 
t{6) = T + B{6 -smO), (7) 

where 9 varies in the range of G [0, 2tt] and T is a constant. For 6 = n, 
the structure reaches to the maximum radius of Vmax = 2^4. Substituting in 
equation ([2]) at the maximum radius provides A = GM/2E which results in 

. 1 Rj 
A = 

6 6i 

On the other hand using equation ([3]) yields: 

A^ = GMB^ 

Substituting equations (EI) and (C!) in ([T]) results the evolution of the density 
contrast in terms of 6 as: 

'^-2n-cos^^)3-l- (8) 



^ Note that the effect of density contrast in the velocity of the structure and sub- 
seque 



sequently on the Kinetic Energy of over-dense region is missed in the text book 
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For the initial condition, considering 5i <^ 1, the initial phase is 9i = 2^/5^. 
Taking T = 0, t{6) will coincide with the cosmic time. From the initial condi- 
tion, B obtain as: 



From equation ([8]) at 6' = 27r/3 the structure enters to the non-linear regime 
(i.e. 6 ^ 1). 

On the other hand for 6 = 27t we have singularity, however before this stage, 
global radial velocity of the structure is converted to the dispersion velocity 
and prevent it from the catastrophic collapsing. The virial theorem provides us 
a radius that the stable condition of the structure is fulfilled. In the next part 
we will extend the top-hat model to the two component fluid of the baryonic 
and dark matter in which they are gravitationally coupled during the evolution 
of the structure. 



2.2 two component top-hat model 



Evolution of the large scale structures indicates that dark matter is an essential 
element for the formation of the structures in the universe. In the standard 
scenario of the structure formation, structures composed by baryonic matter 
and dark matter are in mutual gravitational interaction during their evolution. 
In this section we obtain the evolution of each component in the structure, 
using the top-hat model. 

Let us start with baryonic component. Baryons before the decoupling were 
tightly coupled to the photons and diffusion of the radiation prevents them 
to collapse under their own gravity. The corresponding diffusive mass of the 
baryonic structure with A < l^iff is called the Silk mass with the mass of 
Ms = 6.2 X 10i2(^/Q^)3/2(^/^2)-3/4j^^^ Howcver, after decoupling of 

the baryons from the cosmic microwave background radiation, Jeans length 
decreased rapidly and baryonic structure could grow freely (see Fig. [1]). The 
corresponding Jeans mass after the decoupling is 

kT 

Mj = p,{- f ^ lO^Me, (9) 

GphTJlp 



where all the parameters are calculated at the last scattering surface with 
T ~ SOOO-ft' and nip is the mass of the proton. The rest of the scenario is the 
gravitational interaction of the baryonic matter larger than the corresponding 
Jeans mass with the gravitational potentials that have already been made by 
the dark matter structures. 
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Comparing the mass of galaxies and cluster of galaxies with the Silk mass 
at the last scattering surface shows that we can set Sf, < 10~^ which can be 
ignored compare to the density contrast of dark matter {6h ^ 5^ ~ 10"'^). 

We consider two spherical regions with radii of Rb{t) and Rd{t), and total 
masses of Mj, and for the baryonic and dark matter components of the 
structure. These two spheres are coincided on each other at the initial time, 
but due to the different initial conditions they will evolve with different rates. 
Similar to the first part of this section, the momentum and energy conservation 
equations for the dark matter is written as 



ii..-G'Ji±^, (10) 

We let the initial density-contrast for the baryonic matter {Sb — 0) which 
provides a zero initial peculiar velocity for the baryonic component. Comparing 
to the dark matter sphere, the baryonic sphere will expand faster. According 
to this picture from the dynamics, the dark matter will interact gravitationally 
only with a fraction of the baryonic matter inside the dark matter sphere. 

Similarly, the dynamics of baryonic sphere is given by: 



R,^-G^, (12) 



E,= '-Rl-G^'^ + ^'\ (13) 

2 ^ Rb ^ ' 

As the baryonic sphere is always larger than the dark matter sphere, this 
component interacts with all the dark sphere. The energy of the two spheres 
can be obtained from Eq.([S]). Setting zero peculiar velocity for the baryonic 
sphere, results in Eb = —kf'^dl'^^ and for the dark matter component = 

To calculate the initial condition for the density contrast of dark matter 5\'^'' , 



we use the Harrison-Zeldovich power-law spectrum (jl2l : Il3l ) as 

P{k) = Ak"", (14) 



where we adapt n = 1. The corresponding mass variance of this spectrum is 
cr^ = 4^fc^. An important specification of this spectrum is that at the entering 
time of the structure to the horizon the density contrast has an invariant value 
of aenter = 2TcA^/'^/{9tl). Usiug A ~ (28.6/1"^ Mpcy, the density contrast for 
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the entering time is Senter — 6 x 10~^. Our aim at this stage is applying the 
Harrison-Zeldovich spectrum to have the density contrast of the dark matter 
at the decouphng time. The density contrast at the last scattering depends 
on the epoch that structure enters to the horizon. We divide the evolution of 
density contrast into two area of radiation and matter dominant epochs. In 
these two phases, the dynamics of dark matter changes with different rates in 
terms of scale factor. 

Let us first consider the structures that enter the horizon at the radiation 
dominant epoch. The structures having smaller than the horizon mass at the 
equality time (M < M//(eg)) will enter horizon at the following redshift: 

, M , 1 ^ ^ 

ZenteriM) ^ Z^gi ) ^. (15) 

^^■iH[Zeq) 



Letting horizon mass at the equality epoch Mui^Zeq) ~ 5 x lO^^(f2/i^)~^M0, the 
corresponding entering redshift of a structure to the horizon can be calculated. 
For instance the dark matter component of a structure with a galaxy mass 
enters the horizon at Zeneter{,go.laxy) ~ 5.90 x 10''. The structures grow during 
Zenter (enter time) to Zeq (equality time) by a logarithmic factor of 

8{aeq) ~ 5 ln(^^) 5 (enter), (16) 

O'enter 



where we use Zi,q = 3454T392 (j3|). Structure with a galaxy size at the present 
time will grow by the factor of 20 during z^^nter to z^g period which results a 
density contrast of about 6eq{galaxy) ~ 8.50 x 10^'^ at equality time. Once the 
universe enters to the matter dominate era, structure start to grow propor- 
tional to the scale factor and at the decoupling time the density contrast of 
dark matter will reach to 6{zdec) = idee I Center x Senter- For a structure with the 
galaxy mass, the dark matter density contrastgrows up to 5dec — 2.7 x 10^^ at 
the decoupling time ( we use Zdee = 1088^2 (jsl)). Fig. ([2]) shows the dynamics 
of the radii of dark matter and baryonic components of a galaxy in terms 
of redshift. We use ACDM model for the background with the parameters 
of fi^^ = 0.3, Vt^^ = 0.7 and Hq ~ 70 Km/sMpc~^ . The time dependence 
of density contrast also is shown in Fig. ([3]). The evolution of each compo- 
nent is calculated up to the virilization time. We note that the dark matter 
component virialize earlier than the baryonic part. 

We calculate the evolution of the gravitationally coupled two component in 
the top-hat model, for the structures in the range of 10® to 10^^ solar masses. 
These structures, depending on their masses, will enter the horizon at different 
epochs, smaller sooner and the larger later. Means that the smaller structures 
grow faster than larger ones (see Tableland Fig.H]). This result is in agreement 
with the down-top scenario of the structure formation. 
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An important observational data for examining this model is comparing the 
power spectrum of the large scale structures at the present time with that of 
model. We do this comparison by calculating the mass variance for a sphere 
with the radius of i?, in terms of power spectrum. The total mass Mji{r) in a 
sphere with the radius R, centered on the point r is: 



\r-r'\<R 

= J p{r')e{R - |r - r'\)d\', (17) 

where is the step function. Using the Fourier transformation of density 
contrast and step function, the mass variance is related to the power-spectrum 
through 

/AMr\2 r, . .odk 



M 



J \W{kR)\'\A,\'^, (18) 



where = A;^/^ 1 5^ | /V^vr and W{kR) is the Fourier transformation of step 
function. The function W{kR) cuts off the integral f|T8l) for A; > 1/R and since 
Afc is an increasing function of k, so the integral is generally dominated at 
k ^ 1/R and we have (jl3 ) : 

(^)~A, k^ 1.38/R. (19) 



Figure ([5]) shows the calculated power spectrum in two component top-hat 
model in terms of k. This power spectrum derived from this model has almost 



the same trend as the observations (jl4j ). We should mention that since the 



power-spectrum of galaxies has scale dependent biasing, these data cannot be 
used for comparing with the dark matter spectrum. To have a comparison of 
this model with the observed data, we use erg derived from the gravitational 



weak-lensing observations (llStl.We point out that ag derived from the weak- 
lensing deosn't suffer from the biasing problem and it probes directly the 
distribution of the dark matter. Weak lensing observations provide 0.62 < 
(jg < 1.32 which is compatible with that of top-hat model, erg — 0.8. 

In the next section we will discuss about the effect of background dynamics 
on the evolution of the structures namely the effect of a variable dark energy 
model. 
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3 Top-Hat Model in Variable Dark Energy Background 



In this section we discuss the effect of background dynamics on the growth 
rate of the structures. Recent observations of SNIa and CMB show that uni- 
verse is mainly made by an exotic fluid so-called dark energy which speeds 



up its expansion (Il6l ). The standard solution for interpretation of the posi- 
tive acceleration of the universe is including the cosmological constant to the 
Einstein equation. The best fit with the observations provides fi^-* = 0.3 and 
^a'' =0.7, where f^A is the density parameter corresponds to the cosmologi- 
cal constant. While ACDM model provides a good fit to the SNIa and CMB 
data, however it suffers the coincidence problem and finite tuning of the cos- 
mological constant at the early universe. One of the solutions is considering 
a variable dark energy model. In this model an evolving scalar field generates 
the energy and the pressure of the dark energy and for the later times in the 
history of universe, it provides a positive accelerating universe. 

In this section we study the effect of variable dark energy on the dynamics 
of the spherical collapse model as a function of the redshift. The dark energy 
can influence on the growth of large-scale cosmological structures through (i) 
the background effect, in which the dark energy changes the expansion rate 
of the background and (ii) dark energy may deviate from the homogenous 
distribution due to the gravitational interaction with the dark matter. The 
feedback of the dark energy is changing the growth rate of the dark matter. 



These effects have been studied in the following works (llTl : Il8l : Il9l : l20l : l2ll : |22 



23|) 



In the case of ACDM universe the spherical collapse is similar to that in 
CDM model, except the cosmological constant that changes the growth of the 
structure though altering the background dynamics. For the CDM model the 
virialization radius (Rvir) is the half of the maximum radius of the structure 
(Rmax)- In the ACDM universe the virialization radius is smaller than that of 
in CDM model |l7'). Mota and Van de Bruck (MB) considered the spherical 



collapse for different potentials of the quintessence models (l23l). They showed 
that the predictions of the spherical collapse depend on the dark energy model. 
In this scenario, during the collapse of over-dense regions, the dark matter 
enters the highly non-linear regime while the perturbation in dark energy 
deviates slightly from the background. They also showed that if the dark 
energy equation of state is assumed to be constant, the differences between the 
homogenous and inhomogeneous cases are small. The advantage of considering 
the dark energy in the spherical collapse model is that it predicts that the 
cosmic structures such as the clusters of galaxy have been collapsed prior to 
the epoch of z ~ 1.4, compatible with the observations of the most distant 



cluster (l29t). 
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In this part we examine the effect of a parameterized dark energy on the 
evolution of the structures through the dynamics of the background. We take 
the equation of state from Wetterich (2004): 



where luq is the equation of state at the present time and b is the bending 
parameter (j24j ). The best fit from the comparison of the cosmological data 



with the model resuhs in 6 = LSSl^go and wq = -lAbtoH (1251). The density 



parameter of this dark energy from the continuity equation changes as: 



\3[l+ui(z;6,a;o)] 



(21) 



where is the energy density of dark energy at the present time and 
uj{z]b,ujQ) = uJo/[l + 61n(l + z)] is the average of the equation of state in 
the logarithmic scale. Using the Hubble parameter for the fiat universe, 

H'{z; 6,^o) = H',p^\l + zf + + zf + + zr'^% (22) 



we obtain the dynamics of the scale factor for the various bending parameters 
as shown in Figure ([6]). To calculate the evolution of the baryonic and dark 
matter components of the structure, we use equations ffTOj) and f|T2|) to obtain 
the radius of the structure as a function of time. The evolution of the Hubble 
parameter as a function of redshift is given by equation (l22l) . Eliminating 
time in favor of the redshift for the dynamics of the structures we obtain 
R = R{z) as shown in Fig. Here we fix the value of wq and calculate 
the dynamics of the structures for various bending parameters. Increasing 
the bending parameter causes the structure forms at earlier epochs but with 
smaller radius. Also the effect of Wq on the formation of structures, while b is 
a fixed value, is shown in Fig. ([8]). Decreasing ujq causes the structures form 
faster with smaller radius. 

In the dark energy models, the evolution of the structures deviate from that in 
ACDM model. This effect results from the change in the background dynamics 
due to the dark energy. In addition to this effect, the total energy of the 
structure will be changed due to the dark energy effect. The effect of a variable 
dark energy in contrast to ACDM model is that we can have a none-zero 
contribution of the dark energy at the early epoches of the universe. So for 
a given Hq at the present time, we expect to have larger Hi at the early 
times. Looking to the total energy of a structure E = —3/2HfRf6i shows 
that having larger Hi causes more bounded structure. This effect provides a 
negative energy for the structure which will produce smaller Rmax and R^r as 
well as earlier virialization time to the structure. 
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In the next section we will study the effect of cooling of baryonic matter at 
the finial stage of virialization on structure formation. 



4 Cooling Mechanism 



Once the baryonic structure reaches to the virial condition, radial velocity of 
the structure converts to the dispersion velocity or in anther word the baryonic 
gas thermalizes. The temperature of structure, using the kinetic energy of 
baryonic particles can rise up to lO^K and from the Saha equation we will 
have an ionized medium. The ionized plasma then cools down and the result 
is more contraction of the baryonic structure. On the other hand through the 
cooling, the baryonic structure can fragment into the small parts to generate 
stellar systems. During the cooling of baryonic matter, it will lose its kinetic 
energy and hence falls into the gravitational potential well and subsequently 



gains the kinetic energy (1261 ). This cooling and heating continues until the 



system reaches to a quasi-steady state. 

A simple parameter that represents the cooling of a gas is the ratio of cooling 
to the free fall time scale, r = tcooi/tgrav The cooling time is defined by tcooi = 
E/E 3pkBT/2fiA{T) and the dynamical time also results from the time 
scale for the free fall of a structure and is given by tdyn = |(2GM_R~^)~^/^. 
If r > 1 then the gas can cool; but as it cools the gas can retain the pressure 
support by adjusting its pressure distribution. If r < 1 the gas will cool rapidly 
to a minimum temperature. This will lead to the loss of pressure support and 
the gas will undergo an almost free-fall collapse. In this case fragmentation 
into stellar units can occur. There are various physical processes contribute in 
cooling, depending on the temperature of the plasma. We assume cooling is 
dominated by plasma Bremsstrahlung radiation and recombination. At tem- 
perature larger than lO'^K, Bremsstrahlung dominates the cooling whereas in 
the range of 10*^ ~ 10^, recombination of gas is the main source of cooling. 
The net cooling rate is 



where A(T) is the radiative cooling function and is expressed as 

A{T) = {AbT-^ + AnT^)p' (24) 



Ab oc e^/m^/^ represents cooling due to the bremsstrahlung and Aji — e^nipAB 
arises from the recombination. This expression is valid for temperatures above 
10"^ K. For lower temperatures, the cooling rate drops drastically because H 
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can no longer be significantly ionized by collisions. The radiative cooling func- 
tion for temperatures above 10^ is well approximated by 

A(T) = 2.5 X lO-'^T^-^. (25) 



Here we study the cooling condition of baryonic part of structure after the 
thermalization which r = tcooi/tdyn <^ 1 is hold and we have almost free fall 
condition for the structure. Fig. ([9]) shows the variation of r in terms of redshift 
for a galaxy mass structure, starting from the thermalization time to the quasi 
steady state. The duration of cooling for this structure in the ACDM model 
is about Az ~ 0.07. Fig. flTUl) shows the evolution of the size of structure 
in terms of redhsift for this period. After thermalization structure freely falls 
until it reaches to a quasi steady state. Fig. (fTTl) also indicates the evolution 
of the density contrast including the cooling effect after thermalization time. 

To see the effect of variable dark energy on cooling time, we take a small struc- 
ture with the mass of lO^M©. This mass is suitable for studying the evolution 
of the globular clusters and corresponds to the baryonic Jeans mass after the 



decoupling (l27l ). According to the structure formation scenario we expect that 
this structure has a dark matter component. However the observations of dis- 
persion velocity show that these structures almost have no dark matter. The 
pre-galactic model for the globular cluster and tidal striping of dark halo by 
Galaxy can explain the lack of dark matter in these structures ( jisi ) . It should 
be noted that in the formation of these structures we don't consider merging 
effect as naturally is taken account in the N-body simulations. 

Variable dark energy causes the structure thermalizes and cools at the higher 
redshifts (see table |2]). For instance, taking the parameters wo = —1.45 and 
b = 1.35, structures with globular cluster mass thermalize at 2; ~ 2.17 and stop 
cooling at 2; ~ 1.94. For a larger structure with a galaxy mass of IO^^Mq, the 
corresponding thermalization redshift occurs at z ~ 0.7 (see Fig. [7j). Compar- 
ing with the globular cluster we can conclude that the first star bursts would 
happened in the globular clusters. The high metalicity with the low rotation 
of globular cluster in the galactic halo supports this hypothesis. 



5 Conclusion 



In this work we extend simple top-hat model into two component dark matter- 
baryonic structure. The initial condition for the sub-horizon size structures 
is taken at the last scattering surface. The scale of X < dn for the struc- 
tures guarantees applying approximately the Newtonian mechanics. Using the 
Harrison-Zeldovich spectrum for the perturbations of dark matter, we obtained 
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the evolution of each component of the structure. 

For the dark matter part, we showed that density contrast in the small mass 
structures grows faster than the larger ones and subsequently reaches to the 
maximum radius and virializes at the higher redshifts. This behavior of struc- 
tures implies that the star burst should take place at the smaller structures as 
dwarf galaxy and globular clusters. An observable parameter of structures to 
compare with the model is the power spectrum. We calculated the mass vari- 
ance of structures in various scales at the present time and compared as ~ 0.8 
from the model with that of observation from the weak lensing. For studying 
the dynamical effect of background on the evolution of the structures, we ap- 
plied a logarithmic variable dark energy model and showed that the structures 
in this model evolve faster than that of ACDM. 
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Fig. 1. Size of horizon of universe (short dashed line), size of structure (long dashed 
line) and Jeans length (solid line) of baryonic structures in terms of scale factor in 
logarithmic scale. At the last scattering surface the Jeans mass of baryonic matter 
decreases and baryonic structure with sub-horizon scale can grow after this epoch. 
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Fig. 2. Dependence of radii of halo (dotted-line) and baryonic (dashed-line) compo- 
nents for astructue with a galaxy mass in top-hat model in terms of redshift, com- 
pared with that of background (solid- line) . Background is taken ACDM model with 
the corresponding parameters of = 0.3, Jlj^^ = 0.7 and Hq ~ 70 Km/ sMpc~^ . 




Fig. 3. Density contrast evolution of baryonic (dashed-line) and dark matter (solid- 
line) in terms of redshift. The horizon line represents 5 = 1, separate the liner and 
non- linear regimes. 
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Fig. 4. Dependence of Mass (normalized to M//)to the characteristic redshifts of 
the structures. Sohd hne represents redshift corresponds to the maximum radius of 
a structure in terms of mass. Dashed hne represents the dependence of viriahzation 
redshift to the mass of dark matter structure. Dotted and dashed-dotted hues rep- 
resent the maximum radius and viriahzed redshifts for the baryonic component of 
the structure, respectively. 
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Fig. 5. The power spectrum calculated by the two component top-hat model is 
derived in terms of k. 
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Fig. 6. Dynamics of scale factor in terms of cosmic time (normalized to to) with a 
variable dark energy model given by Eq. ()20|) . Solid line stands for 6 = (ACDM), 
dashed line b = 0.5, dotted line b = 1.0 and dashed-dotted line b = 1.5. 
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Fig. 7. The effect of bending parameter in logarithmic variable dark energy model 
on the dynamics of structure in terms of redshift. Red lines represent the baryonic 
component and blue lines stand for the dark matter component. The bending pa- 
rameters are chosen as 6 = (solid line), b = 0.5 (dashed line), b = 1 (dotted line) 
and b = 1.5 (dashed-dotted line). The equation of state is fixed with loq = —1.45. 
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Fig. 8. The effect of equation of state in logarithmic variable dark energy model 
on the dynamics of structure in terms of redshift. Red lines represent the baryonic 
component and blue lines stand for the dark matter component. The equation of 
states are chosen as wq = (solid line), wq = —0.5 (dashed line), wq = —1 (dotted 
line) and wq = —1.5 (dashed-dotted line). The bending parameter is fixed with 
b= 1.35. 
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Fig. 9. Variation of r in term of redshift for a galaxy mass structure, r is plotted 
from thermalization till the quasi-steady state phase of baryonic structure. Dashed 
era represents the stable zone for the structure where the cooling is negligible. 
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Fig. 10. Dynamics of radius of galaxy mass structure in terms of redshift (dashed 
line). After thermalization we will have free fall collapse of structure until quasi 
steady state time. Solid line represents the dynamics of background for comparison. 
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Fig. 11. Evolution of density contrast in terms of redshift for baryonic (dashed 
line) and dark matter (solid line) components of a galaxy mass structure. After 
thermalization we will have almost four order of magnitude increase of the baryonic 
density contrast. 
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Table 1 

Numerical results from the evolution of structures with various masses. First column 
shows the range of mass of structures and subscript in M represents mass in 10"Mq . 
Second column is the corresponding redshift to entering to the Horizon. The density 
contrast of dark matter at the decoupling epoch is in third column. Forth and Fifth 
columns are the corresponding redshift of maximum radius for the dark matter and 
virialization of the dark structure. Sixth and seventh columns are the same as the 
forth and fifth columns for the baryonic structure. 
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Table 2 

The effect of bending parameter on duration of cooling of baryonic structure in log- 
arithmic variable dark energy model. Mass of structure is taken 10^ Mq. At the first 
column we fixed {wq = —1.45) and the second column shows the duration of cooling 
for various bending parameters. In the third column we fixed bending parameter 
(6 = 1.35) while the equation of state 'u;o has different values with corresponding 
duration of redshift is indicated in the fourth column. 
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-1.45 
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[2.12, 1.90] 


b = 


1.50 


[2.25,2.00] 


LVo = 


-1.50 


[2.19,1.95] 
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